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ABSTRACT 

Context. Doppler tomography is a technique that has revolutionised the interpretation of the phase-resolved spectroscopic observations 
of interacting binary systems. 

Aims. We present the results of our investigation of reversing the velocity axis to create an inside-out Doppler coordinate framework 
with the intent to expose overly compacted and enhance washed out emission details in the standard Doppler framework. 

Methods. The inside-out tomogram is constructed independently of the standard tomogram by directly projecting phase-resolved 
spectra onto an inside-out velocity coordinate frame. For the inside-out framework, the zero-velocity origin is transposed to the outer 
circumference and the maximum velocities to the origin of the velocity space. We test the technique on a simulated system and two 
real systems with easily identifiable features, namely the accretion disc and bright spot in WZ Sge, and spiral shocks in IP Peg. 
Results. Our tests show that there is a redistribution of the relative brightness of emission components throughout the tomograms, 
i.e., where the standard framework tends to concentrate and enhance lower velocity features towards the origin, the inside-out velocity 
framework tends to concentrate and enhance higher velocity features towards the origin. Conversely, the standard framework disperses 
and smears the higher velocities farther away from the origin whereas the inside-out framework disperses and smears the lower 
velocities. In addition, the projection of the accretion disc in velocity space now appears correctly orientated with the inner edge close 
to the maximum velocity origin and its outer edge closer to the zero velocity outer circumference. Furthermore, the gas stream and 
secondary star are projected on the outside of the disc with the bright spot of the stream-disc impact region on the disc’s outer edge in 
the inside-out velocity space. 

Conclusions. We conclude that inside-out Doppler tomography complements the already powerful tomographic techniques in the 
analysis of spectroscopic emissions from interacting binary systems. 
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1. Introduction 

Cataclysmic variables (CVs) are interacting binary systems in 
which a white dwarf (primary) accretes material from a late-type 
main-sequence star (secondary). The secondary fills its Roche 
lobe and material flows through the inner Lagrangian point (Li) 
towards the primary. The material may form an accretion disc 
around the primary before it is finally accreted (see Warner 1995, 
for a comprehensive review of CVs). 

Marsh & Horne (1988) developed Doppler tomography as a 
technique aimed at constructing a two-dimensional velocity im¬ 
age (tomogram) of the accretion discs of CVs using an emis¬ 
sion line in its spectra sampled at a number of orbital phases. 
Doppler tomography has revolutionised the interpretation of or- 
bitally phase-resolved spectroscopic observations of interacting 
binary systems. It has become a valuable tool for resolving the 
distribution of line emission in CVs and other binary systems 
(see, e.g.. Marsh 2001, Astrotomography workshop review). Sig¬ 
nificant contributions have been, amongst others, the discovery 
of two-arm spiral shocks in the accretion disc of the dwarf nova 
IP Peg (Steeghs et al. 1997) and mapping of the accretion stream 
in the polar HU Aqr detailing the transition from ballistic to mag¬ 
netic flow (Schwope et al. 1997). Recent literature also abounds 
with examples where emission from the secondary and the bright 


spot of the stream-disc impact region in CVs have been resolved 
in Doppler tomograms, such as CTCV J1300-3052 (Savoury 
et al. 2012) and V455 And (Bloemen et al. 2013). 

We present a complementary technique for an inside-out ve¬ 
locity projection for Doppler tomography. This paper expands 
on the preliminary results presented in Kotze & Potter (2015). In 
Section 2 we first review the cartesian coordinate frame used in 
Doppler tomography and introduce a polar coordinate frame for 
the new inside-out velocity space. Next, in Section 3 we compare 
the standard and inside-out Doppler tomograms for a simulated 
disc-accreting system and two real CV systems, WZ Sge and IP 
Peg. Also in Section 3, we investigate how the inside-out frame¬ 
work redistributes relative contrast levels amongst the various 
emitting components of observations of a variety of CV features. 
Finally, in Section 4 we provide a summary and conclusions. 

2. Velocity space 

2.1. Velocity in cartesian coordinates 

We first review the cartesian coordinate frame that co-rotates 
with the binary system, as described by Marsh & Horne (1988). 
The top panel of Fig. 1 shows the cartesian spatial coordinates 
(normalised by the binary separation a) for a model CV with 
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Fig. 1. Standard spatial and velocity cartesian coordinates for a model 
CV. The C.O.M. of the binary is marked with a plus (+) while that of 
the primary and secondary are marked with crosses (x). Overlays of 
the Roche lobe of the primary (dashed line), the accretion stream (solid 
line) from the secondary towards the primary (up to an azimuth angle of 
30°), as well as radii at decrements of 0.25 times the outer disc radius are 
shown (red dashed circles). The outer disc radius is the 3:1 resonance 
radius while the inner disc radius is based on an absolute radial velocity 
of 1800 km s“^ There are also two arbitrary reference points, A and B, 
shown in the disc. 


an accretion disc. The assumed binary system parameters for 
the model CV are: a primary mass of O.SSM©, a mass ratio of 
q = 0 . 2 , an orbital period of 0.05669 d (82 min) and an inclina¬ 
tion of 11°. The origin is placed at the centre of mass (C.O.M.) of 
the binary (for a low q system this is very close to the C.O.M. of 
the primary), the A-axis runs along a line through the centres of 
the primary and secondary, the F-axis runs along a line paral¬ 
lel to the velocity vector of the secondary and the Z-axis along a 
line through the C.O.M. of the binary perpendicular to the orbital 
plane. In CVs where Roche-lobe overflow is considered to be re¬ 
sponsible for the accretion stream and disc, these features are pri¬ 
marily conflned to the orbital plane, therefore a two-dimensional 
layout in the orbital plane, i.e., the AF-plane, is sufficient. The 
assumption that all motion is parallel to the orbital plane forms 
one of the postulates of Doppler tomography (Marsh 2001). The 
binary orbital phase zero is defined as the point of mid-eclipse 
of the primary by the secondary. The orbital motion is counter¬ 
clockwise around the C.O.M. of the binary. 

Every spatial point (x, y) in the system has a two- 
dimensional cartesian velocity coordinate (Vj^, V 3 ;) as shown in 
the bottom panel of Fig. 1. Given the definition of the spa¬ 
tial frame it follows that, for example, the C.O.M. of the pri¬ 
mary is moving in a negative F-direction (V 3 ; is negative), but 



Fig. 2. Velocity bins in the cartesian velocity framework (i.e., in the 
lower right quadrant of the Doppler map in Fig. 1). The figure shows 
the positions of a lower and a higher velocity pixel (red and blue, re¬ 
spectively). 


with no movement in the A-direction (Vx is zero). Similarly, the 
C.O.M. of the secondary is moving in a positive F-direction (v^ 
is positive), but also with no movement in the A-direction (Vx 
is zero). Structures such as the accretion stream and disc are 
more complex since they have a rotational velocity about the 
C.O.M. of the binary and a velocity due to structural motion. 
The increasing velocity of the accretion stream along its ballistic 
trajectory is clearly seen in the absolute increase in its compo¬ 
nent. Due to the Keplerian nature of the accretion disc its veloc¬ 
ity profile appears reversed, i.e., the higher velocities of the inner 
edge of the disc are farther from the origin in velocity space and 
conversely for the lower velocities of the outer edge of the disc. 
The Keplerian nature of the disc is further emphasised by the 
radii at decrements of 0.25 times the outer disc radius. The sec¬ 
ondary and ballistic stream also appear to be ‘inside’ the disc in 
velocity coordinates. 

We note that it is not possible to invert the velocity map to 
produce a spatial map, i.e., every velocity coordinate does not 
have a unique spatial coordinate making the inversion from ve¬ 
locity to spatial coordinates impossible without additional infor¬ 
mation (e.g., Harlaftis et al. 2004). In reality, orbitally phase- 
resolved spectroscopic observations are mapped directly onto 
the velocity map. Effectively each pixel, centred on a velocity 
coordinate (Vj^, ^y) in the velocity map, represents a velocity bin 
with a width (dvx, dvy). In the cartesian coordinate frame of 
square pixels, dvx = dVy for all velocity bins throughout the 
map. For example, in Fig. 2 we show a red and a blue pixel 
in the lower right quadrant of the Doppler map in Fig. 1 with 
velocities (v^^, ^y) = (150, -150) and (690, -1710) kms“\ re¬ 
spectively. Both pixels have the same velocity bin widths of 
dvx = dVy = 60kms“^ corresponding to a range in velocity 
magnitude of 85kms“^ between opposite velocity bin corners. 
However, the range in velocity angle {dO) is 22.6° and 2.4°, re¬ 
spectively, i.e., the range in velocity angle decreases as a function 
of radial distance from the origin. Consequently, the lower radial 
velocity bins (pixels) sample a larger velocity vector range than 
the higher radial velocity bins (pixels). 

The sinusoidal radial velocity curve ( 0 ), traced by each ve¬ 
locity coordinate (Vj^, V 3 ;), is a function of the orbital phase 0 and 
centred on the mean (systemic) velocity y of the binary, which 
is given by: 

Vr ( 0 ) = y - Vjc cos 27r0 -f V 3 ; sin Incp. ( 1 ) 
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Fig. 3. As for Eig. 1 but in polar coordinates. 

To simplify the discussion this can be equally given in polar co¬ 
ordinates as: 

Vr (0) = y -V cos(6 27r0), (2) 

where v is the velocity magnitude and 6 is the velocity angle. 
The range in velocity magnitude of 85 km s“^ between opposite 
pixel (velocity bin) corners results in four different sinusoidal 
amplitudes. The range in amplitude is the same for both the 
red and blue pixels. However, the lower velocity red pixel sub¬ 
tends a larger velocity angle (6) which corresponds to a larger 
range in the relative phases of the four sinusoids and therefore 
bound a larger area in the radial velocity-orbital phase space. 
Conversely, the blue pixel has a narrower bounding area in the 
radial velocity-orbital phase space. The implications of this ef¬ 
fect is discussed in later sections. 

2.2. Velocity in polar coordinates 

We continue our discussion in the polar coordinate system as this 
is more conducive to the circularly symmetric nature of Doppler 
tomograms and gives a more direct indication of velocities and 
directions. The same CV model shown in Fig. 1 is shown in 
Fig. 3, but in spatial and velocity polar coordinates. Spatial co¬ 
ordinates (normalised by the binary separation a) are now given 
as the radial distance r from the C.O.M. of the binary and the 
polar angle 6s (not to be confused with 6 in the velocity frame) 
measured in an anti-clockwise direction from the line between 
the C.O.M. of the binary and the secondary (Fig. 3 top panel). 
The corresponding velocity map in polar coordinates is shown 
in the bottom panel of Fig. 3, again with the same model CV as 
for Fig. 1. Velocity magnitude v increases as a linear function of 


90 ° 




Fig. 4. Inside-out velocity polar coordinates. The top and bottom panels 
show the model CV centred on the anti-C.O.M. of the binary and the 
anti-C.O.M. of the primary, respectively. See text for full description. 

distance from the origin with a direction given by the polar angle 
6 in an anti-clockwise direction measured from the line drawn 
from the C.O.M. of the binary horizontally to the right. All the 
CV model parameters and the arbitrary points A and B are the 
same as in the previous section. Therefore, for example, the sec¬ 
ondary has a velocity magnitude of 400kms“^ with a velocity 
direction (angle) of ^ = 90° (i.e., = 0 and V 3 ; = 400kms“^ in 

cartesian coordinates). 

2.3. Velocity in inside-out polar coordinates 

Having described the velocity polar coordinate frame in the pre¬ 
vious section, we now investigate constructing an inside-out ve¬ 
locity polar coordinate frame with zero velocity on the outer cir¬ 
cumference and the maximum velocities around the centre of the 
coordinate frame. The velocity magnitude v now increases as a 
linear function of distance from the zero velocity outer circum¬ 
ference towards the origin with a direction given by the polar 
angle 6. In order to make the inside-out framework complemen¬ 
tary to the standard framework, and to facilitate a direct compar¬ 
ison between the frameworks, the default origin of the inside- 
out framework is set to the maximum radial velocity used to 
construct the standard framework. Effectively this is the maxi¬ 
mum radial velocity extracted from the phase-resolved spectra 
of a source emission (or absorption) line, i.e., the ‘edge’ of the 
data projected onto the framework to construct the tomogram. 

The inside-out velocity map of the model CV presented in 
the previous figures is shown in the top panel of Fig. 4. Essen¬ 
tially the zero velocity origin (corresponding to the C.O.M. of 
the binary) is now transposed to the outer circumference and 
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Fig. 5. Velocity bins in the inside-out polar velocity framework (i.e., in 
the lower right quadrant of the Doppler map in the top panel of Fig. 4). 
The figure shows the positions of a lower and a higher velocity pixel (red 
and blue, respectively). These pixels are the inside-out polar equivalent 
of the cartesian pixels shown in Fig. 2. 


the maximum outer velocities to the central origin, giving the 
new inside-out velocity framework. The central origin, however, 
represents a discontinuity and the inner 4 velocity bins (pixels) 
around it are therefore ignored. We would like to point out that 
the spatial map (Fig. 3 top panel) remains unchanged. As one 
can see (Fig. 4 top panel) the inner and outer edges of the accre¬ 
tion disc are now correctly orientated. Also the ballistic stream 
is now on the outside of the disc and ‘curves’ inwards as it ac¬ 
celerates towards the disc and primary. The Roche lobe of the 
secondary is also now on the outside of the disc albeit upside 
down. It is upside down because it is orbiting as a solid body 
with the outside moving faster than the inside, as opposed to the 
Keplerian velocities of an accretion disc and the ballistic veloc¬ 
ities of the accretion stream. The Roche lobe of the primary is 
now the outer bounded circular ring area, i.e., between zero ve¬ 
locity and the dashed line, that also contains the C.O.M. of the 
primary (indicated by a x). 

The Keplerian velocity radial profile of the disc is still vis¬ 
ibly apparent as can be seen by the equidistant disc radii con¬ 
tours (Fig. 4 top panel). We explored replacing the linear velocity 
axes with a Keplerian velocity (oc 1 / ^/r) profile in an attempt to 
produce a more equidistant velocity contour profile. However, a 
correct calculation of the Keplerian profile would require knowl¬ 
edge of the primary mass of the actual system under investiga¬ 
tion. This is rarely known and would introduce a new parameter 
that would no longer give a map that is simply a direct projection 
of the observations. 

We note that the inner edge of the accretion disc is not cen¬ 
tred on the origin (Fig. 4 top panel). This is because the disc is 
centred on the C.O.M. of the primary, whereas the origin repre¬ 
sents the anti-C.O.M. of the binary, i.e., the C.O.M. of the bi¬ 
nary is effectively the outer circumference of zero velocity with 
its anti-C.O.M. at the origin. Hence the Keplerian velocities of 
the disc will be offset by an amount corresponding to the orbital 
motion of the primary. This offset is also present in standard ve¬ 
locity maps but is more apparent in the inside-out framework. 
The bottom panel of Fig. 4 shows the same model but with the 
orbital motion of the primary subtracted from the model velocity 
profile. This effectively places the primary on the zero velocity 
outer circumference and its anti-C.O.M. at the origin. The effects 
of centring Doppler tomograms on the C.O.M. of the binary or 
of the primary are investigated in the next section. 


As described in Sect. 2.1 each velocity bin (square pixel) 
has the same velocity magnitude range; however, the velocity 
vector range decreases as a function of radial distance from the 
central origin. To demonstrate this for the inside-out framework 
we show in Fig. 5 the equivalent red and blue pixels shown in 
Fig. 2. In the inside-out framework the higher velocity pixels 
are closer to the origin and therefore sample a larger velocity 
vector range than in the standard framework and conversely for 
the lower velocity pixels. This can be seen clearly by comparing 
Figs 2 and 5. The implications of this effect are discussed in the 
following sections. 

3. Doppler tomography: standard and inside-out 
projections 

For real CVs, the orbital and structural velocities can be ob¬ 
served in the Doppler-shifted variations in the multi-component 
emission (and absorption) lines of their optical spectra. The or- 
bitally phase-resolved spectroscopy of individual spectral lines 
can be projected onto a predefined velocity framework to pro¬ 
duce a Doppler tomogram. The interpretation of such tomograms 
would then proceed by comparison with overlays of model ve¬ 
locity profiles of the likes described in the previous section (e.g.. 
Marsh et al. 1994; Wolf et al. 1998). Traditionally, Doppler to¬ 
mograms are produced by a projection onto a cartesian coordi¬ 
nate frame, but given the circularly symmetric velocity profile 
of binary systems we elected to use a polar coordinate frame for 
our standard and inside-out tomograms as discussed in earlier 
sections. 

We have modified the fast maximum entropy Doppler to¬ 
mography code presented by Spruit (1998) to incorporate the 
inside-out velocity framework^. With the modified code the ob¬ 
served Doppler-shifted spectra are projected directly onto an 
inside-out velocity coordinate frame. All the standard and inside- 
out Doppler tomograms presented hereafter have been created 
using the unmodified and modified code, respectively. 

3.1. Simulated system: accretion disc, bright spot and 
secondary 

The results of applying standard and inside-out Doppler tomog¬ 
raphy to a simulated data set are shown in Fig. 6. A single line 
emission spectrum was constructed to simulate a disc-accreting 
system with an orbital period of 0.05669 d (82 min), a mass ra¬ 
tio of ^ = 0.2, a primary mass of 0.85Mo and an inclination of 
77°. The emission line was made of four components: two rel¬ 
atively narrow Gaussian profiles to represent emission from the 
accretion disc’s bright spot and the irradiated face of the sec¬ 
ondary star, a broad double peaked profile to represent emission 
from an accretion disc and some artificial noise. Appropriate 
phase-dependent Doppler shifts were then calculated for each 
of the components before summing to produce the trailed spec¬ 
tra shown in the top left panel of Fig. 6. The effect of inclination 
was also taken into account. 

The binary system parameters listed above were used to cal¬ 
culate the model velocity profile overlay shown in all the tomo¬ 
grams. Included in the overlay are the Roche lobes of the pri¬ 
mary (dashed line), the secondary (solid line) and a single parti¬ 
cle ballistic trajectory (solid line) from the Li point up to 40° in 
azimuth around the primary, as well as the 3:1 resonance radius 
as the outer edge of the accretion disc (large solid line circle). 
Also shown is an inner Keplerian radius (dot-dashed line) based 

^ Our code may be obtained upon request. 
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Fig. 6. Doppler tomography of synthetic spectra. The standard tomogram and two inside-out tomograms (with the anti-C.O.M. of the binary and 
of the primary at the origin, respectively) are shown from left to right in the middle panels for comparison (see text for full description). The top 
and bottom panels show the input and reconstructed trailed spectra, respectively, for the corresponding tomogram. 


on an absolute radial velocity of 1800kms“^ This represents 
the maximum absolute radial velocity seen in the emission line 
variations. 

3.1.1. Standard Doppler tomogram 

The middle left panel of Fig. 6 shows the standard Doppler 
tomogram. The ring-like feature corresponds to emission from 
the disc with additional emission regions at the locations of the 
bright spot and secondary star. The lower left panel is the recon¬ 
structed trailed spectra from the Doppler code. 

3.1.2. Inside-out Doppler tomogram 

Anti-C.O.M. of binary at origin: The same trailed spectra (top 
middle panel) were used to produce the inside-out Doppler to¬ 
mogram (middle panel). As expected from our modelling in pre¬ 
vious sections, emission from the secondary star is now located 
on the outside of the disc and the bright spot on the disc’s outer 
edge. Furthermore the disc has the ‘correct’ orientation with a 
radial profile that has changed from a more even distribution in 
the standard tomogram, to a more inner disc concentration. This 


can be understood in terms of our discussion in Sect. 2, i.e., in 
the standard tomogram the higher velocity disc emission is dis¬ 
tributed over more pixels covering a larger surface area. In the 
inside-out tomogram the higher velocity disc emission is con¬ 
centrated into fewer pixels closer to the origin and hence appears 
brighter. This is perhaps more consistent with what one might 
expect from disc emission, i.e., the emitting volume in the inner 
edges are much smaller that the outer edges of the disc. How¬ 
ever, the lower half of the disc appears to contain more emission 
at higher velocities than the upper half. This asymmetric artefact 
arises because the Keplerian motion of the accretion disc orbits 
the primary which in turn orbits the C.O.M. of the binary, i.e., 
the disc is projected off-centre. This asymmetry is also present in 
the standard tomograms but is enhanced in the inside-out tomo¬ 
grams due to the brightening of the high-velocity pixels around 
the origin. 


Anti-C.O.M. of primary at origin: We next subtracted the orbital 
velocity of the primary Ki ~ 91 kms“^ (e.g.. Marsh et al. 1987) 
from the phase-resolved spectra (top right panel). As expected, 
the double peak emission from the disc has been ‘straightened’ 


Article number, page 5 of 9 































































A&A proofs: manuscript no. ioDopplerAstro-ph 


in the trailed spectra, i.e., with the orbital velocity of the pri¬ 
mary removed, emission from a given radius in the disc shows 
a constant radial velocity in the trailed spectra. In the recalcu¬ 
lated inside-out tomogram (middle right panel) the C.O.M. of 
the primary now effectively becomes the outer circumference of 
zero velocity with its anti-C.O.M. at the origin. Consequently, 
the disc is now centred on the origin and the asymmetry seen in 
the disc emission has disappeared as expected. 

3.2. WZ Sge in quiescence: accretion disc and bright spot 

WZ Sge is a prototypical member of the CV subclass referred 
to as WZ Sge-type dwarf novae (Bailey 1979). Doppler tomog¬ 
raphy of WZ Sge (e.g.. Spruit & Rutten 1998) has revealed a 
prominent accretion disc and bright spot where the ballistic ac¬ 
cretion stream impacts the outer edge of the disc. This makes it 
an ideal test case for inside-out Doppler tomography. 

The standard and inside-out Doppler tomography results for 
WZ Sge in quiescence are shown in Fig. 7 (similar layout to 
Fig. 6). The tomography is based on the H a emission line from 
spectroscopic observations included with the fast maximum en¬ 
tropy Doppler tomography code (Spruit 1998). The model veloc¬ 
ity profile overlay shown in all the tomograms was calculated us¬ 
ing the orbital period of 0.056687846 d (82 min) determined by 
Patterson et al. (1998), the inclination of 77° obtained by Spruit 
& Rutten (1998) and the mass parameters preferred by Steeghs 
et al. (2007), i.e., a primary mass of 0.85Mo and a mass ratio 
q = 0.09. The overlay includes the Roche lobes of the primary 
(dashed line), the secondary (solid line) and a single particle bal¬ 
listic trajectory (solid line) from the Li point up to 45° in azimuth 
around the primary, as well as the 3:1 resonance radius as the 
outer edge of the accretion disc (large solid line circle). An inner 
Keplerian radius (dot-dashed line) based on an absolute radial 
velocity of 1600 km s“^ which represents the maximum absolute 
radial velocity seen in the emission line variations, is also shown. 

3.2.1. Standard Doppler tomogram 

The most distinct features in the standard tomogram shown in 
the middle left panel of Fig. 7 are the relatively bright accre¬ 
tion disc and the bright spot of the impact region between the 
ballistic stream and disc. The disc appears brighter around the 
assumed outer edge (velocities lower than 1000kms“^) while 
it becomes more diffused towards the assumed inner edge (ve¬ 
locities greater than 1000 km s“^). The extended nature of the 
bright spot (stream impact region) in velocity space is consid¬ 
ered to be the result of the mixing of stream and disc mate¬ 
rial with different velocities (Marsh et al. 1990). The model 
stream lines up well with the brightest part of the impact region 
located at (820kms“\ 150°). There is only a slight enhance¬ 
ment in the emission at the expected velocity of the secondary 
(400kms“\ 90°). 

3.2.2. Inside-out Doppler tomogram 

Anti-C.O.M. of binary at origin: In the inside-out tomo¬ 
gram (middle panel) the disc is now correctly orientated with 
brighter disc emission occurring at generally higher velocities 
(> 1000km s“^). Due to the effect of the ofiT-centre projection 
of the disc there is a slight enhancement in the brightness of the 
disc in the lower half. Although given WZ Sge’s low mass ratio 
this effect is fairly minimal. Similar to the standard tomogram 
the bright spot (stream impact region) also dominates the bright¬ 


ness scale of the inside-out tomogram. However, it is now lo¬ 
cated towards the disc’s outer edge in velocity space and has ac¬ 
quired an extended enhancement in the lower left quadrant cov¬ 
ering (1000km s“\ 180 - 250°). The brightest part of the im¬ 
pact region it is now located at (1050kms“\ 165°), a slightly 
different velocity position than in the standard tomogram. This 
shift is caused by the redistribution of the relative contrast lev¬ 
els throughout the tomogram. The model stream, however, still 
lines up well with the brightest part of the impact region. The 
secondary is not discernible within the low-velocity emission 
covering (0-700 km s“\ 0 - 360°). This is expected since it was 
barely discernible in the standard tomogram and is projected into 
more pixels in the inside-out tomogram. However, by subtracting 
the axisymmetric average of the emission at each radius around 
the origin (not shown), there is a patch of diffuse emission at 
the expected location of the secondary. Similar to the simulated 
system in the previous section, the high-velocity (2000 km s“^) 
noise present at the edge of the spectra is enhanced by the bright¬ 
ening of the high-velocity pixels around the origin. 

Anti-C.O.M. of primary at origin: The velocity of the primary 
Ki ~ 47kms“^ with a phase zero offset 0.12 (Steeghs et al. 
2007) was subtracted to obtain the phase-resolved spectra in the 
top right panel of Fig. 7. Given the relative low velocity of the 
primary that was subtracted, the ‘straightening’ of the double 
peak emission from the disc in the trailed spectra is less pro¬ 
nounced compared to the case of the simulated system in the 
previous section. In the recalculated inside-out tomogram (mid¬ 
dle right panel) the anti-C.O.M. of the primary is now at the 
origin. With the effect of the off-centre projection of the disc re¬ 
moved a more uniform circularly symmetric appearance of the 
disc is recovered. The extended enhancement of the bright spot 
in the lower left quadrant covering (1000kms“\ 180 - 250°) is 
also less pronounced. 

3.3. IP Peg in outburst: spiral shocks 

IP Peg is an eclipsing member of the dwarf nova subclass of 
CVs. Outbursts lasting 10-12 days occur every 70 - 100 days 
during which it brightens by 2 - 3 magnitudes. The spiral shocks 
observed in the accretion disc of IP Peg during an outburst 
(Steeghs et al. 1997) favour the disc instability model (Osaki 
1974) as the cause of the recurrent outbursts. Since the disc ex¬ 
pands during the outburst the outer disc will experience an in¬ 
creased gravitational attraction from the secondary. This causes 
perturbations in the circular Keplerian orbits of the disc material 
in the outer disc, leading to the formation of spiral arms in the 
disc. The spiral structure in the accretion disc of IP Peg makes it 
an ideal test case for inside-out Doppler tomography. 

Standard and inside-out Doppler tomography results based 
on the Hen (44686 A) emission line from the 1996 November 
spectroscopic observations taken during outburst maximum, are 
shown in Fig. 8 (similar layout to Fig. 7 and Fig. 6). These ob¬ 
servations were presented first by Harlaftis et al. (1999) and the 
reader is referred to that paper for a detailed description of the 
data set and reduction procedures. All the tomograms in Fig. 8 
show only the non-axisymmetric emission. The binary system 
parameters as calculated by Copperwheat et al. (2010), i.e., an 
orbital period of 0.1582061029 d (228 min), an inclination of 
83.8°, a primary mass of 1.16Mo and a mass ratio q = 0.48 were 
used to calculate the model velocity profile overlay shown in all 
the tomograms. Included in the overlay are the Roche lobes of 
the primary (dashed line), the secondary (solid line) and a single 
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Fig. 7. Doppler tomography of WZ Sge in quiescence. Same layout as Eig. 6. 
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particle ballistic trajectory (solid line) from the Li point up to 
45° in azimuth around the primary, as well as the tidal radius as 
the outer edge of the accretion disc (large solid line circle). Also 
indicated is an inner Keplerian radius (dot-dashed line) based 
on an absolute radial velocity of 1200kms“\ representing the 
maximum absolute radial velocity seen in the emission line vari¬ 
ations. 

3.3.1. Standard Doppler tomogram 

In the standard tomogram shown in the middle left panel of 
Fig. 8 the emission from the leading side of the irradiated sec¬ 
ondary is seen as a bright spot at (200kms“\ 95°) which ap¬ 
pears to be ‘inside’ the disc. There is also a low-velocity compo¬ 
nent visible at (0 - 300kms“\ 140 - 250°). This low-velocity 
component has been seen in other dwarf novae during outburst, 
but its origin is not understood (Steeghs et al. 1996). The ori¬ 
gin of the high-velocity component seen as a diffuse patch of 
emission at (1000 - 1500kms"\ 0°) is also uncertain. The to¬ 
mogram is dominated by the emission from the shocks in the ac¬ 
cretion disc that creates an extended two-armed spiral structure. 
The two spiral shocks have a clear asymmetry in both brightness 
and velocity. The observed velocity asymmetry of the material 
in the shocks suggests a highly non-Keplerian flow in the disc 


since circular rings of emission are expected from a purely Ke¬ 
plerian disc. The shocks appear to be spiralling ‘outwards’ from 
lower to higher velocities. The upper spiral shock extends to ve¬ 
locities lower than our velocity profile overlay of the outer disc 
radius. This supports the suggestion by Steeghs et al. (1997) that 
the disc possibly expanded beyond its tidal radius early on in the 
outburst as a result of it being triggered by a disc instability. 

3.3.2. Inside-out Doppler tomogram 

Anti-C.O.M. of binary at origin: IP Peg’s mass ratio (0.48) is 
fairly high. Consequently the C.O.M. of the primary is signif¬ 
icantly offset from the C.O.M. of the binary system. This can 
be clearly seen in the offset of the model overlay w.r.t. the ori¬ 
gin in the inside-out tomogram (middle panel). As expected 
the secondary is seen as an extended diffuse patch of emission 
‘outside’ the disc. However, the low-velocity component visi¬ 
ble at (0 - 300kms“\ 140 - 250°) in the standard tomogram 
is no longer discernible. The high-velocity patch of emission at 
(1000 - 1500kms“\ 0°) has become more enhanced and com¬ 
pact given its new projection closer to the origin. The two-armed 
spiral structure associated with the shocks in the disc is still 
clearly present in the inside-out tomogram. However, the asym¬ 
metry in the brightness of the shocks caused by the off-centre 
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Fig. 8. Doppler tomography of IP Peg in outburst. Same layout as Fig. 7 and Fig. 6. Only the non-axisymmetric emission is shown in the 
tomograms. 


projection of the disc is more pronounced than in the standard 
tomogram. The upper shock is more extended and diffuse while 
the lower shock is more compact and brighter. 


Anti-C.O.M. of primary at origin: The phase-resolved spec¬ 
tra in the top right panel of Fig. 8 were obtained by subtract¬ 
ing the orbital velocity of the primary Ki ~ 152kms“^ (aver¬ 
age from Copperwheat et al. 2010). Given the complex struc¬ 
ture of the trailed spectra (Harlaftis et al. 1999) the ‘straight¬ 
ening’ of the double peak emission from the disc is primarily 
seen in the blue peak. These ‘straightened’ spectra were used as 
input to produce the tomogram in the middle right panel. The 
model overlay clearly displays the removal of the C.O.M. off¬ 
set. The secondary has become more compact and brighter. This 
is because the position of the secondary has moved signifi¬ 
cantly closer to the origin leading to its enhancement. Similarly 
the patch of high-velocity emission at (1000 - 1500kms“\ 0°) 
has become slightly more enhanced. Although the origin of 
this feature is unknown it is associated with the low bright¬ 
ness emission seen in the input trailed spectra at high radial 
velocity (> 1000km s“^) between phases 0.4 and 0.8. In ad¬ 
dition an extended low-velocity component is now visible at 


(0 - 250kms“\ 110 - 150°). This low-velocity component is 
part of the low-velocity component of unknown origin (Steeghs 
et al. 1996) covering (0 - 300kms“\ 140 - 250°) in the stan¬ 
dard tomogram. The asymmetry between the brightness of the 
two shocks caused by the off-centre projection of the disc has 
also been removed. The spiral nature of the shocks is clearly ev¬ 
ident and is more intuitive to interpret in respect that it shows 
how they begin at the outer edges of the disc and curve inwards 
towards the primary as they increase in velocity. 


4. Summary 

We have investigated the use of an inside-out velocity coordinate 
frame for Doppler tomograms of non-magnetic cataclysmic vari¬ 
ables. In the inside-out framework zero velocity is on the outer 
circumference and the maximum velocities are around the centre 
of the coordinate frame. We modified the fast maximum entropy 
Doppler tomography code presented by Spruit (1998) to incor¬ 
porate the inside-out framework which allowed us to construct 
inside-out tomograms independently of standard tomograms by 
directly projecting the observed spectra onto the inside-out ve¬ 
locity coordinate frame. 


Article number, page 8 of 9 

































































E. J. Kotze et al.: Exploring inside-out Doppler tomography 


We applied our new inside-out velocity projection to real 
data of the dwarf nova system WZ Sge during a normal faint 
state. The accretion disc has the correct orientation with the 
bright spot located towards the disc’s outer edge (in velocity 
space). In addition there is a redistribution of disc brightness, 
i.e., in the standard tomogram the disc appears to brighten for 
velocities lower than 1000 km s“^ whereas with the inside-out 
velocity projection the disc brightens for velocities higher than 
1000 km s“^ The difference arises because the line emission 
from the higher velocities emanating from the inner regions of 
the disc are spread out over a larger area with more pixels in 
the standard tomogram, whereas the higher velocity emission is 
concentrated in a smaller area with fewer pixels in the inside- 
out tomogram as is more in line with a real disc brightness pro¬ 
file. Given the relative low mass ratio {q ~ 0.1) of WZ Sge this 
real data set only marginally highlighted the vertical asymmetry 
in the brightness distribution in tomograms with an off-centre 
projection of the accretion disc, i.e., with the C.O.M. or anti- 
C.O.M. of the binary at the origin. The asymmetry is somewhat 
enhanced in the inside-out velocity framework (anti-C.O.M. of 
the binary at the origin), but it can be removed by subtracting 
the orbitally induced Doppler velocity of the primary from the 
phase-resolved input spectra. 

We also applied our new projection technique to real data of 
the dwarf nova system IP Peg during an outburst which shows 
emission from the irradiated face of the secondary and a dis¬ 
tinct spiral structure in the accretion disc. The secondary is now 
placed ‘outside’ the correctly orientated disc in velocity space. 
The dominant spiral structure of the shocks in the disc is well 
preserved and the shocks are now spiralling ‘inwards’ from 
lower to higher velocities. Due to the relative high mass ratio 
~ 0.5) of IP Peg the vertical asymmetry in the brightness dis¬ 
tribution in the tomograms with the C.O.M. or anti-C.O.M. of 
the binary at the origin is more pronounced than in the case of 
WZ Sge. This implies that it is crucial to consider the mass ratio 
{q) of the system when interpreting the brightness distribution in 
a tomogram. The asymmetry in the brightness distribution is re¬ 
moved in tomograms with the anti-C.O.M. of the primary at the 
origin resulting in a more circularly symmetric brightness distri¬ 
bution. This, together with the ‘inwards’ spiralling appearance 
of the shocks, creates a velocity image which is perhaps more 
intuitive with the expected radial disc profile of the shocks. 

We conclude that our new technique of inside-out Doppler 
tomography is complementary to the existing technique. The 
standard projection used in Doppler tomography tends to con¬ 
centrate and enhance lower velocity features while higher veloc¬ 
ity features are more separated and dispersed. Conversely, the 
inside-out velocity projection tends to concentrate and enhance 
higher velocity features while lower velocity features are more 
separated and dispersed. This is perhaps more consistent with 
the emission distribution in disc-accreting CVs where the high- 
and low-velocity emissions are produced in smaller and larger 
fractions of the total system volume, respectively. 
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